Stark effect on the exciton spectra of vertically coupled quantum dots: horizontal field 

orientation and non-aligned dots 
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We study the effect of an electric-field on an electron-hole pair in an asymmetric system of 
vertically coupled self-assembled quantum dots taking into account their non-perfect alignment. We 
show that the non-perfect alignment does not qualitatively influence the exciton Stark effect for the 
electric field applied in the growth direction, but can be detected by application of a perpendicular 
electric field. We demonstrate that the direction of the shift between the axes of non-aligned dots 
can be detected by rotation of a weak electric field within the plane of confinement. Already for a 
nearly perfect alignment the two-lowest energy bright exciton states possess antilocked extrema as 
function of the orientation angle of the horizontal field which appear when the field is parallel to 
the direction of the shift between the dot centers. 



PACS numbers: 73.21. La,71.35-y,71.35.Pq,73. 21. Fg 



I. INTRODUCTION 

The effect of the electric field on the exciton recombina- 
tion energy is used to probe the properties of the electron 
and hole confinement in separate self-assembled quantum 
dotsii^ through a strong deformation of the carrier wave 
functions. For vertically coupled quantum dots^ the elec- 
tric field oriented in the growth direction leads to a redis- 
tribution of the carriers between the dots.— '^'^'^ Typical 
electric field applied for vertically coupled quantum dots 
is of the order of 20-30 kV/cm, i.e., by an order of magni- 
tude smaller than the one used to probe the confinement 
potential^ of a single dot. Therefore, such experiments 
on the clcctric-ficld effect for coupled dots resolve rather 
the properties of the molecular coupling than the fabri- 
cation dependent details of the confinement in separate 
dots. Moreover, the charge redistribution induced by the 
electric field between the different dots have naturally 
a much stronger energy effect than deformation of the 
wave functions inside each of the dots. Due to these fea- 
tures, previous simple modeling® (done in parallel with 
the experimental work'*) using a square well vertical and 
harmonic lateral confinement^ led to the correct descrip- 
tion of the observed*' - electric field induced dissociation 
of the electron-hole pair in an avoided crossing of dark 
and bright energy levels It also successfully predicted- 
non-trivial features of the photoluminescence spectrum 
observed during the negative trion dissociation by the 
electric field. ^ 

A vertical electric field was used very recently^ to verify 
the growth process of intentionally strongly asymmetric 
double quantum dots. In this paper we demonstrate that 
the horizonal electric field (perpendicular to the growth 
direction) can be useful to estimate the non-perfect verti- 
cal alignment of the dots and to determine the direction 
of the horizontal shift between them. To account for the 
horizontal field effects we replaced the harmonic lateral 
profile of our previous model® by a quantum well and de- 



veloped a computational approach for treating excitons 
in systems with no axial symmetry and for an arbitrary 
electric field direction. 

The vertical electric field-induced redistribution of 
charge carriers is smooth only for the electron which 
tunnels much more effectively between the dots and is 
responsible for the most characteristic experimental fea- 
tures of the spectrum4ii^i^iiS Similarly, it is the electron 
redistribution between the dots due to a rotation of a 
weak horizontal field that allows for the detection of the 
non-perfect alignment, the hole charge being only shifted 
within the separate dots. Therefore, we decided to use 
for simplicity the single band model to describe the hole. 

The paper is organized as follows. Section II describes 
the computational approach, our numerical results are 
given in Section III, and Section IV contains the summary 
and conclusions. 



II. DOUBLE DOT MODEL AND METHOD OF 
CALCULATIONS 



We work in the single valence band approximation, for 
which the Hamiltonian of the electron-hole pair is 
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with Te, Yh - electron and hole coordinates, rrie, - ef- 
fective band masses for the electron and the hole, e the 
dielectric constant, and F is the electric field. The elec- 
tron confinement potential of the double dot system is 
taken in the form of a sum of double disk-shaped quan- 
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FIG. 1: The electron confinement potential [Eq. (2)] in meV, 
in the y — plane, for the ofi'set between the axes of the 
dots d — 10 nm (equal to the dot radius R) and the vertical 
distance between the dot centers of h — 10 nm. The height 
of both the dots is 2Z = 4 nm. The lower (negative z) dot 
has depth V;*^ = 518 meV, the dot for positive z has depth of 
V,f = 508 meV. 



turn wells with smooth walls (see Fig. [T]) 



Veir) 



(2) 



-K7(i + ( 
-K7(i + 



{x-d/2y+y' 



)")(!+ (^)") 



here R stands for the radius of a single dot, Z is half of 
its height, the dots centers are placed symmetrically with 
respect to the origin at points {x,y,z) = ±{d/2,0,h/2) 
(the spacer layer is then equal to 6 = ft. — 2Z), Vj and 
Vg" are the lower and upper dot depths, respectively. For 
simplicity we take the same size of both dots, the in- 
evitable asymmetry between them (due to size difference, 
indium composition, nonsymmetric strain, etc.) can for 
the limited purpose of the present paper be effectively 
taken into account by applying different depths of the 
dots. The confinement potential for the hole is of the 
same form, only with depths V^, replacing and 
. The bottom of the conduction band of the barrier 
material (GaAs) is taken as the reference energy for elec- 
trons, and the top of the GaAs valence band as the refer- 
ence level for the holes, i.e., the eigenvalues of (1) have to 
be shifted up by the GaAs band gap to give the photon 
energy measured in a luminescence experiment. 

We previously^ demonstrated, that for an interacting 
electron-hole pair both the energy spectrum and the par- 
ticle distribution are nearly the same when the asymme- 
try of the double well is exclusively introduced for one of 
the particles, since it is translated by the interaction into 
an asymmetry for the other particle. In this paper we 
assume that the dot situated on the negative side of the 
z = plane is deeper by 10 meV for both the electron 



and the hole = ^ 10 meV, = + 10 meV, and 
= 508 meV, = 218 meV is taken the same as in 
our previous paper— assuming that the dots embedded in 
a pure GaAs matrix are made of Ga2;Ini_a;As alloy with 
X — 0.66. Accordingly^ we take me — 0.037, m/i = 0.45, 
and e ~ 12.5. The size parameters are taken from Ref<^, 
diameter 2i? = 20 nm and height 2Z = 4 nm (therefore 
b = h — A nm). 

The Schroedinger equation with Hamiltonian (1) is 
solved using the exact diagonalization (configuration in- 
teraction) approach keeping a complete account of the 
electron-hole correlation. As the basis for the electron- 
hole pair we take products of single-electron states and 
single-hole states 



■jk 



(3) 



where /i"'', fj^^ are the single-electron and single-hole 
wave functions of state j and k, respectively. The single- 
particle wave functions are obtained in the multicenter 
basis of displaced Gaussian functions, 

M 

fir) = ^c,exp[-ai{{x-x^)'^ + {y-yif) 

i 

-Pr[z~Z,f], (4) 

where Ci are the linear variational parameters, ai and [3i 
are the nonlinear variational parameters describing the 
strength oii—th Gaussian around point {xi, yi, Zi), whose 
coordinates are variational parameters by themselves. 
We took 11 basis functions for each of the dots: one 
situated at the dot center ±(c?/2, 0, /i/2), eight around 
it within the plane of confinement Zi = ±/i/2, on a 
circle with optimized diameter. Two additional cen- 
ters are placed above and below the center of each dot 
±(d/2, 0,h/2 + ±A) with A optimized variationally. Ad- 
dition of more centers outside the plane of confinement 
or including more centers within it does not improve the 
results significantly. In the configuration basis (3) we in- 
clude all the single-particle wave functions obtained in a 
diagonalization of the electron and the hole Hamiltonians 
each in a basis of form (4) (with position of centers and 
other nonlinear parameters optimized separately for the 
electron and hole), which eventually yields 484 localized 
basis functions to treat the exciton. The recombination 
probability for the exciton state with wave function cf) is 
calculated from the envelope wave function as 



P 



d^r(f>ire,rh)S^{re-Th) 



(5) 



Two-dimensional version of the multicenter basis was 
previously used for laterally coupled dotsi^i Our studyi^ 
on electron systems in a single circular dot performed 
with a similar technique showed that the superposition 
of Gaussians very well approximates the angular momen- 
tum eigenstates. 
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FIG. 2: The exciton spectra for the electric field oriented along the z direction for a vertical distance between the dots centers 
h — 10 nm for (a) perfectly aligned dots d = 0, (b) horizontal distance between the vertical symmetry axis of the dots d = 5 
ran, and d — 10 nm (c). 
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FIG. 3: The electron and hole probability density distribu- 
tions for the ground state and the first excited bright state in 
the absence of the electric field for d = 5 nm and h = 10 nm 
(a), and h = 8 nm (b). 



III. RESULTS 

The effect of the unperfect alignment of the dots for 
the exciton Stark effect due to the electric field oriented 
in the growth direction is presented in Fig. [21 The fig- 
ure shows the dependence of the exciton spectra on the 
vertically oriented electric field for the barrier thickness 
b = 6 nm (ft, = 10 nm). For the dots aligned vertically 

(a) as well as for the offset between the axes of d 5 nm 

(b) and d = 10 nm (c) at F = we observe two bright 
lines (thickness of lines is set proportional to the recom- 
bination probability). The particle distribution for the 
bright states is displayed in Fig. [31 Fig. [3^a) shows the 



cross section {y = 0) of the electron and hole probability 
density distributions (obtained through integration of the 
two-particle probability density in the coordinates of the 
other particle) for both the bright states at F = 0. We 
notice that in these states the hole is entirely localized in 
one of the dots, and the electron tends to accompany the 
hole. In the ground-state the electron is entirely localized 
in the deeper (lower) dot. In the excited bright state, the 
electron is predominantly localized in the upper dot, but 
a small part of the electron wave function leaks to the 
other (deeper) dot. For thinner spacer layer b = A nm 
{h = 8 nm) a leakage of the electron density to the dot 
without the hole appears in the ground state [Fig. [Sfb)] , 
and in the excited state it is strengthened. Note, that in 
the excited state the hole density maximum is displaced 
to the left within the upper dot following the electron 
density noticeably shifted to the lower dot. 

The vertically oriented electric field F > tends to 
push the electron to the upper (shallower) dot. In the 
ground state near F = 20 kV/cm we see [cf. Fig. [HJa-c)] 
an avoided crossing of the bright energy level with both 
carriers in the deeper dot with a dark energy level of sep- 
arated carriers (hole stays in the lower dot, the electron 
is transferred to the upper one). The crossing is avoided 
due to the electron tunnel coupling between the dots. In 
the upper bright state - no avoided crossing is observed 
for F > 0, the positive electric field stabilizes the elec- 
tron in the upper dot. On the other hand, for the upper 
bright energy level, the negative electric field of smaller 
absolute value (around 10 kV/cm) is enough to transfer 
the electron from the upper to the lower (deeper) dot 
which results in an avoided level crossing similar to the 
one appearing in the ground-state [see the upper avoided 
crossing at left of F = in Fig. [HJa-c)]. 

We see that for the displaced axes of the dots [Figs. [H 
(b-c)] the dependence of the spectra on the vertical elec- 
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FIG. 4: The exciton spectra for the electric field oriented along the x direction for the vertical distance between the dots 
centers h = 10 nm for (a) perfectly aligned dots d = 0, (b) horizontal distance between the vertical symmetry axis of the dots 
d = 5 nm, and (c) d = 10 nm. Single-particle densities for the bright energy levels labeled by numbers 1-4 in (b) are presented 
in Fig. 6. 




FIG. 5: Electron and hole densities for 3 lowest-energy states 
for d=0 nm and ft = 10 nm taken within the plane of con- 
finement of the lower, deeper dot {z — —5 nm). For the field 
of _F = 20 kV/cm parallel to the x axis (arrows at the elec- 
tron and hole panels at the top-left side of the figure show the 
direction of the electric force acting on the electron and the 
hole). The energy levels of these states are displayed in Fig. 
BTa). 



trie field is qualitatively similar to the perfect alignment 
case [Fig. ^a)] but with an avoided level crossings width 
which is smaller due to the suppressed electron tunnel 
coupling as the barrier thickness in increased. 

Pronounced qualitative differences of the spectra for 
aligned and not aligned dots appear in case the elec- 
tric field is oriented horizontally. Fig. SJa) shows the 
spectrum for perfectly aligned dots. For both the bright 
states (with both carriers localized in the same dot) the 
electric field separates the electron and the hole within 
each dot pushing the carriers at its opposite sides, which 
leads to a decreased recombination probability and a de- 



creased energy. The first excited state, which is dark and 
twice degenerated at F = 0, corresponds to both carriers 
in the deeper dot, but with a hole excitation (of p type). 
The horizontal electric field lifts the degeneracy of the p 
energy level of the hole. The electron (left panel) and the 
hole (right panel) density for F = 20 kV/cm applied in 
the X direction is shown in Fig. [5]for the ground-state and 
the two lowest-energy excited states with p-excitations of 
the hole. The plot was made in x,y coordinates for the 
plane of confinement of the lower dot z = —h/2 [all the 
three states correspond to carriers totally localized in the 
lower dot cf. upper panels of Fig. [2Ia)]. The lower dark 
energy state with p-hole excitation has a nodal surface at 
the plane y = 0, and the parity with respect to this plane 
is conserved when the field is applied in the x-direction 
(see the panel of Fig. [5] for the first excited state), hence 
the zero recombination probability for all F. The other p- 
level, higher in energy becomes bright at larger F, when 
the hole parity with respect to the direction perpendic- 
ular to the field is destroyed (see the panel of Fig. [5] for 
the second excited state). We note that in all the three 
states, for F = 20 kV/cm the electron density is not per- 
turbed by the horizontally applied field, but its effect on 
the hole (for which the confinement is weaker) is clearly 
visible. 

For non-perfectly aligned dots the component of the 
electric field parallel to the direction of their relative hor- 
izontal displacement redistributes the carriers between 
the dots. The low-energy spectrum plot for d = 5 nm 
{h still equal to 10 nm) contains 4 bright energy levels 
denoted by 1, 2, 3, and 4 in Fig. Hfb). The modification 
of the carrier distribution by the electric field in these 
states is shown in Fig. [51 In the ground state [for F = 
both carriers in the lower dot - see Fig. [3fa)] the elec- 
tron passes to the right (shallower) dot near F — +60 
kV/cm. On the other hand, in the excited state which is 
bright at F = and marked by '2' in Fig. [5] both carri- 
ers are localized in the upper shallower dot. The electron 
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FIG. 6: The y = cross section of the electron and hole probability density distributions for the ground state and the first 
excited bright state for different values of the electric field in the a;-direction (given at the top of the figure for each set of plots) 
for d = 5 and ^ = 10 nm. The selected bright states denoted by numbers 1-4, as in Fig. \MP)- The contour plots for states 1 
and 2 at F = were shown in Fig. |3] 




FIG. 7: The exciton spectra for the electric field oriented 
along the x direction for a vertical distance between the dots 
centers h — 8 nm for (a) horizontal distance between the 
vertical symmetry axis of the dots d = 5 nm. (c) and d = 10 
nm. 



FIG. 8: The ground-state mean electron and hole position 
dependence on the horizontal electric field oriented along the 
X direction for h = 8 nm (a) perfectly aligned dots and (b) 
offset of the dots axes of d = 5 nm. The energy levels for 
parameters used in (b) are displayed in Fig. [3 a). 



is transferred to the left dot (lower and deeper) for the 
field around F — —30 kV/cm. Both these electron trans- 
fers (ground state at -t-60 kV/cm and excited bright state 
at -30 kV/cm) are associated with avoided crossings of 
bright and dark energy levels. These avoided crossings 
become smaller for an offset of 10 nm between the axes 
of the dots, the case presented in Fig. [HJc). 

Fig. [7] shows the spectra for stronger interdot cou- 
pling (h = 8 nm). Compared to /i = 10 nm case, the 



weak-field extrema of the energy levels are shifted more 
distinctly off the F — Q point due to a larger value of 
the built-in dipole moment resulting from the shift of the 
electron density related to the tunnel coupling [see Fig. 
[3fb) and notice opposite shifts for the extrema of two 
lowest bright energy levels and much larger value of the 
dipole moment for the excited state] . The dipole moment 
(/i = |e| (< r/i > — < Te >) is proportional to the differ- 
ence in the mean positions of the electron and the hole. 
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FIG. 9: The exciton spectra for the electric field F = 20 kV/cm oriented in the x — y plane as function of the angle between 
the field and the x axis (F = F(cos((/)), sin((/)), 0)). The plots in the upper (lower) row correspond to /i = 10 nm {h = 8 nm). 
Figures (a-b), (c-d), (e-f) and (g-h) correspond to d = 2, 5, 10 and 15 nm, respectively. Electron and hole densities for energy 
levels of plot (d) are displayed in Fig. 1101 




FIG. 10: The y = cross section of the electron (the panels 
marked by 'e') and hole (the panels marked by 'h') probability 
density distributions for the ground state (denoted by 1 in a 
square frame) and the first excited bright state (denoted by 
2) for d = 5 and h = 8 nm. The corresponding energy levels 
are displayed and marked by the same numbers in Fig. HJd). 
The horizontal axis corresponds to x variable, the vertical - to 
z variable. The electric field is kept at |F|=20 kV/cm, each 
panel corresponds to a different orientation of the electric field 
within the x — y plane given by ij) which is the angle between 
the F and the x axis. The left plots in each column show the 
electron distribution and the right ones-the hole distribution 
The plots for states 1 and 2 at F = were shown in Fig. |3jb) 
(ground and excited state, respectively). 



The impact of the horizontal field (a;) on the spatial po- 
sition of particles in the ground-state is plotted for /i = 8 
nm in Fig. [S] For perfect alignment d = [Fig. [5Ja)] 
the field applied in x direction does not affect the z po- 
sition of the particles, although one can notice a shallow 
minimum for the electron at F = 0. The horizontal field 
slightly strengthens the electron confinement within each 
of the dots thus enhancing the interdot vertical coupling 
and weakly shifting the electron towards the upper dot. 
Fig. [3^ a) shows a stronger reaction of the horizontal po- 
sition of the hole than of the electron, which was already 
noted in the context of Fig. O For non-perfect alignment 
[d = 5 nm, see Fig. MJ^)] the horizontal field (x) leads to 
a very strong dependence of the vertical positions {z) of 
the particles. For F > the electron is transferred to the 
dot at right and is shifted to the top while hole is stabi- 
lized in the left dot. Due to the electron tunnel coupling 
the electron position is continuous and a smooth function 
of F. On the contrary the ground-state values of the hole 
position have a strong jump at F < —40 kV/cm, where 
the crossing of energy levels appear in the ground-state 
[cf. Fig. [7ja)]. Note that also the electron positions are 
modified during this jump - the electron, localized in the 
lower (left) dot tends to follow the hole when it leaves to 
the upper (right) one. 

Fig. [5] shows the low-energy exciton spectrum for a 
fixed, rather small, length of the electric field of F = 20 
kV/cm as function of the angle </> that it forms with the 
X axis, F = F(cos(0), sin((/()), 0). The upper panel of the 
figure corresponds to /i = 10 nm and the lower to h = 8 
nm. The plots from left to right are calculated for off- 
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sets d = 2,5, 10 and 15 nm (for perfectly aligned dots 
d = the spectra are independent of </>) . Already for the 
smallest offset [Figs. [9] (a-b)] we notice that the ground 
state energy is minimal when the field is oriented to the 
right {(p = 0), i.e., when it tends to transfer the elec- 
tron to the shallower right dot, and maximal when the 
field is oriented to the left (tends to keep the electron in 
the deeper - left dot) . For the excited bright state (when 
both carriers tend to remain in the shallower - right dot), 
the energy is maximal when the field pushes the electron 
to the right (reduced penetration of the electron to the 
other - deeper dot) and minimal in the opposite case 
(enhanced tunnehng to the other dot). The energy de- 
pendence of the excited bright state on is stronger than 
for the groimd state. This is due to the larger electron 
tunnel coupling for the excited state [see Fig. [3]. The 
electron is simply more willing to pass to the deeper dot. 

For a larger offset oi d — 5 nm [Fig. Hfc-d)] the depen- 
dencies of the bright energy levels on the angle </> preserve 
their character but become more pronounced. The elec- 
tron and hole densities in the y = plane are plotted in 
Fig. [ini for the two lowest-energy bright energy levels of 
Fig. md) and for opposite field orientation. They are to 
be compared to the F = case displayed in Fig. [HJb), 
which shows that in the absence of the field the proba- 
bility to find the electron in the upper dot for the ground 
state was smaller than the probability to find the elec- 
tron in the lower one in the excited bright state. Field of 
20 kV/cm at (/) = (F oriented in x direction) reverses 
this relation, as the field shifts more of the electron to 
the upper dot when in the ground state and removes the 
electron from the lower dot when in the excited state. A 
shift of the hole distributions to the left and the squeeze 
of the distribution to the left side of the dot is also vis- 
ible. For the electric field oriented along the y direction 
{(j> = 7r/2) the electron distribution of the carriers be- 
tween the dots is similar to the F = case. For the field 
oriented antiparallel to the x axis {4> = tt) the electron 
is almost entirely localized in the left dot when in the 
ground state and almost completely so when in excited 
state. 

The anti-locking of the energy extrema as function of 
(j) for the two lowest bright states observed in Figs. [9] 
(a-d) at the electric field orientation matching the direc- 
tion of the shift between the dots can be explained in 
the following way. The electric field leads to a decrease 
of the energy levels through the separation of the elec- 
tron and hole charges [see Eq. (1)]. In the ground state 
[see Fig. [10] the field at the angle = {+x direc- 



tion) enhances the separation stimulating the electron to 
leave the hole [hence the energy minimum see Fig. 9(d)], 
and at (j) — TT {—x direction) prevents the electron-hole 
separation (hence the maximum). In the excited bright 
state the effect of the field on the carrier separation is 
opposite. For the excited state the modulation of the re- 
combination probability is more pronounced, the electron 
coupling and the electron switching between the dots is 
stronger, hence the more pronounced energy dependence 
on [Fig. 9(d)]. 

For the offset of dot axes equal to the radii d ^ 10 nm 
- Figs. ISUe-f) the excited bright energy level is involved 
in an avoided crossing with a dark energy level. For even 
larger offset d = 15 nm - Figs. [9l(g-h) an avoided crossing 
appears also in the ground state. The appearance of 
these anticrossings is due to the stronger energetic effect 
of constant F for increased d. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have shown that for vertically coupled 
dots their non-perfect alignment does not qualitatively 
influence the exciton Stark effect of the electric field ori- 
ented in the growth direction. Although for perfectly 
aligned dots the horizontal electric field only deforms the 
electron and hole density within each of the dots, for 
non-perfect alignment it leads to a redistribution of the 
particles between the dots, i.e., in the direction perpen- 
dicular to the field. We demonstrated that due to the 
relatively strong electron tunnel coupling and strong elec- 
tron confinement within each of the dots the horizontal 
electric field can be used to tune the electron distribution 
in the vertical direction only weakly affecting the electron 
charge distribution within each of the dots. On the other 
hand, due to a weaker hole confinement the horizontal 
field distinctly affects the hole distribution within each of 
the dots, and the hole transfer between the dots is rapid 
due to negligible tunnel coupling. We have shown that 
a rotation of the electric field vector within the plane of 
confinement can be used to detect the non-perfect align- 
ment. To determine the direction of the horizontal shift 
between the non-aligned dots, and to estimate the length 
of the shift, they can be observed from the exciton pho- 
toluminescence spectrum. 
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